Supersolid phase of strongly correlated bosons in an optical cavity 
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We numerically simulate strongly correlated ultracold bosons in a high-finesse optical cavity by 
means of Bosonic Dynamical Mean Field Theory. The complete phase diagram is established, which 
contains two novel self-organized quantum phases, supersolid and checkerboard solid, in addition to 
conventional phases such as superfluid and Mott insulator. At finite but low temperature, thermal 
fluctuations are found to enhance the buildup of the self-organized phases. We demonstrate that 
cavity-mediated long-range interactions can give rise to stable supersolid and checkerboard solid 
phases even in the regime of strong s-wave scattering. In the presence of a harmonic trap, we 
discuss coexistence of these self-organized phases, as relevant to experiments. 
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Experimental realizations of atomic many-body sys- 
tems coupled to a high-finesse cavity have recently at- 
tracted a large amomit of attention In particular, 
the self-organized phase of atoms induced by coherent 
scattering between pump laser and cavity mode has been 
predicted theoretically [2| , and confirmed experimentally 
for laser-cooled atoms in a transversally pumped cav- 
ity 0. However, only recently it has become possible 
to combine a high-finesse cavity with an ultracold quan- 
tum gas in the strong-coupling regime and to experimen- 
tally investigate properties of a Bose-Einstein conden- 
sate (BEC) in an optical cavity [l-Q- ^ phase transition 
from a normal to self-organized phase in an open system 
has been realized Q, and a lifetime up to 10ms of the 
self-organized phase has been achieved which indicates a 
steady state. Up to now, these experiments have, how- 
ever, focused only on weakly interacting condensates. On 
the theory side, there is a lack of quantitative predictions 
for strongly-correlated bosons coupled to an optical cav- 
ity, even though an extended Bose-Hubbard model has 
been derived [1, [13] which describes the ultracold gas 
trapped in a periodic optical potential generated by the 
high-finesse cavity. Recently, theoretical studies of the 
BEC-cavity system have predicted that the ground state 
can be Mott-insulating with finite photon-excitations of 
the cavity mode [ll|, ll2|. However, the robustness of this 
self-organized phase against strong contact interactions, 
finite temperature and the inhomogeneity induced by an 
external trap remains an important open issue. 

To bridge this gap, here we numerically investigate 
the buildup of self-organized phases in ultracold bosonic 
gases coupled to a single-mode cavity field, pumped by 
a laser beam in the transverse direction. This setup is 
similar to a two-dimensional (2D) classical optical lat- 
tice but with a quantized field in the cavity direction. 
Since the cavity field mediates long-range interactions be- 
tween atoms 3" l3| , we investigate the system by means 
of real-space bosonic dynamical mean field theory (RB- 
DMFT) which captures both strong correlations and spa- 



tial inhomogeneity as well as arbitrary long-range order 
in a unified framework 16|. 

Motivated by the recent experiment Q, we consider a 
system of ultracold ®^Rb atoms with s-wave scattering 
length ds = 5.77 nm and atomic transition wavelength 
A = 780.2 nm, which is driven by a linearly polarized 
standing-wave laser with a red-detuned wave length Xp ~ 
784.5 nm in the direction perpendicular to the cavity 
axis. The setup of our simulation consists of the optical 
cavity in the x-direction, driven by a pump laser in the z- 
direction, and a strong confinement freezing the motional 
degree of freedom of the atoms in the third direction 17 1 . 
We choose the cavity decay rate as k = SOOuir which is 
close to the experimental value of k = 27r x 1.3 MHz 
where cofi is the frequency corresponding to the recoil 
energy, Er = Hlor = h^/{2mXl)i^ 2ii x 3.8 kHz). We 
choose the light shift as Uo = Qo/^a = —0.1, which 
leads to an atom-cavity coupling strength go two orders 
of magnitude larger than the cavity decay rate k and thus 
implies that the system is in the strong-coupling regime 
of cavity QED [7|, where Aa denotes the atom-pump 
detuning. 

This system can be described by an extended Bose- 
Hubbard model [^, [13]. We further assume the cavity 
mode to be in a coherent state to simplify the atom-cavity 
coupling, which is in good agreement with experimental 
results [8|. Within this approximation, the cavity mode 
is described by a complex amplitude a, and the param- 
eters of the extended Bose-Hubbard model only depend 
on the average photon numbers. We thus finally obtain 
the effective Hamiltonian employed in the following cal- 
culations: 
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where bl (bi) denotes the bosonic creation (annihilation) 
operator for a Wannier state at site i. Here (Jz) is 
the effective nearest-neighbor hopping amplitude in the 
X- [z-) direction, with the hopping in a;-direction de- 
termined by the cavity mode, /i is the effective chem- 



ical potential, Vi 



with the strength Vtrap of 



the external harmonic trap, and U = Airash^ /m is the 
Hubbard interaction strength. The cavity mode am- 
plitude a = r,,sJ[,Z^i~^Y{b^b^)/{A'^ + iK) ^ with 
= A, - Uo{J§J:^{^b^) + JfE(,,>(^l^.»,/s deter- 
mined self-consistently by the density distribution of the 
atoms. Jg, Jq and Jf denote the onsite single-particle 
matrix elements of the potential generated by the cav- 
ity mode, by scattering between pump laser and cavity 
mode via single atoms, and the first-order tunneling ma- 
trix element between nearest-neighbor sites of the cavity 
mode standing wave, respectively. 7/off = — ■\/|VJ,C/o| de- 
notes the effective pump strength into the cavity through 
atomic scattering, Ac the cavity-pump detuning, and 
Vp the depth of the standing-wave potential created by 
the pump laser in the z direction. The hopping ampli- 
tudes in X- and z-direction for nearest neighbors are given 
by X,z/Er = (4/V^)(K,./i;fl)(^/^^exp(-2v/y.,J^fl) 
and the Hubbard interaction parameter by U/En = 
4v^(a,/Ap)(Ky,14/i;3)(i/4) ^here K is 

the optical lattice depth in the x- {y-,z-) direction. For 
the onsite coupling matrix elements we use a Gaussian 
approximation of the Wannier states. 

The main challenge now is to determine the steady 
state of the BEC-cavity system described by the Hamil- 
tonian ([1]) in the Wannier basis. Here we apply real- 
space bosonic dynamical mean-field theory (RBDMFT) 
[16| which provides a non-perturbative description of the 
bosonic many-body system both in three and two spa- 
tial dimensions (considered here) [l^. RBDMFT, which 
is capable of including the inhomogeneity of a trapped 
system as well as strong correlations between the atoms, 
assumes the self energy to be local but site-dependent. 
In our calculations, we choose the recoil energy En (w^) 
as the unit of energy, and set h = 1. 

We first investigate the robustness of the super- 
solid phase against interactions for experimentally rel- 
evant parameters. The supersolid is characterized by 
coexistence of the staggered order parameter $ = 
{E^h^yblb,)/{J2,blb^) and superfluid order = (6). 
There are two possible signs of $, i.e. the majority of 
the atoms occupy even sites for <!> > or odd sites for 
$ < 0. Intuitively, if the pump laser is strong enough 
to stabilize a larger atom density at the even sites, and 
at the same time we choose a negative shifted cavity de- 
tuning < 0, this implies that the coherent scattering 
between the pump laser and the cavity mode generates 
a potential with minima at the even sites, as indicated 
by the staggered term of equation ([T|. As a result, the 
corresponding potential will attract more atoms toward 
even sites and the system self-organizes into a steady 



state. In the following, we will confirm this heuristic 
argument via numerical simulations based on RBDMFT. 
To ensure that the tight-binding approximation is valid, 
we add an external optical lattice in the cavity direction 
with a depth of Vcxt = 5Er. 



-500 



a: -1000 
LU 

u 
< 



-1500 - 




-2000 



-2500 

10 12 14 16 18 20 

Vp/ER 

FIG. 1. Phase diagram at zero temperature for a total particle 
number A^tot ~ 506 and two different scattering lengths (fis = 
5.77 nm). The two phases are superfluid and supersolid. The 
cavity decay rate is set to n — SOOEr and the light shift is 
Uo = —O.IEr. Inset: critical strength Vp of the standing- 
wave pump laser vs. temperature at flxed cavity detuning 
A, = -IOOOEr. 

Fig. [T] shows the phase diagram of the system with 
scattering length I.Oos and f.25as at zero temperature 
for a square lattice with TViat = f 6 x f 6. We find two dif- 
ferent phases which are characterized by $. For a weak 
pump laser, the system is superfluid with homogeneous 
density distribution and $ = 0. In this case, the mean 
photon number in the cavity is zero. On the other hand, 
if the pump laser is strong enough, more photons are 
scattered into the cavity mode and the atoms organize 
themselves into a checkerboard pattern with |<I>| > 0. 
Our simulations thus clearly confirm the existence of the 
supersolid phase for single-component Bose gases in the 
cavity in the presence of strong onsite interactions. In 
Fig. [T]we observe that the trend of the phase boundary 
from superfluid to supersolid as a function of Vp at fixed 
scattering length is consistent with experiment 
The phase boundary is considerably shifted upwards for 
larger scattering length, which indicates that more pump 
laser power is needed to drive the system into the self- 
organized phase. We also observe that onsite interactions 
have a more pronounced effect on the buildup of the su- 
persolid phase for a stronger pump laser field. Generally, 
there exists also an unstable state for positive shifted cav- 
ity detuning A^ > [s, 21|, which is beyond the scope of 
this work. 

We also investigate the effect of finite temperature on 
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p at low but finite temper- 
since thermal fluctuations excite the atoms from 



the critical pump strength, as shown in the inset of Fig. [T] 
We observe a minimum of 
ature, 

the ground state and thus reduce the energy gap between 
the homogeneous and the self-organized state. As a re- 
sult, less power of the pump laser is needed to stabilize 
the supersolid. On the other hand, at high temperature, 
thermal fluctuations tend to smear out the self-organized 
density pattern, and as a result, more power is needed 
to stabilize it. Note that the long-range order </> 7^ 
at T > in two dimensions is a mean-field artefact in 
the thermodynamical limit, while in reality, th e sy stem 



should exhibit a Kosterlitz-Thouless transition [22 1 
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FIG. 2. Properties of the self-organized phases of strongly 
interacting bosons in the homogenous lattice. The blue curve 
corresponds to the filling dependence of checkerboard order 
"1? at zero temperature, the green (red) triangle denotes the 
point where the system is in the Mott insulator (checker- 
board solid), (a)-(d): density distribution of superfluid, su- 
persolid, Mott insulator and checkerboard solid, respectively 
in real space (left) and in quasi-momentum space (right), 
corresponding to the densities marked by the red arrows 
in the main figure. Other parameters are Ac — —500Eh, 
K = SOOEr, Uo = -0.1 Er and Vp = 15Er. Inset: Mehing 
of the supersolid phase with increasing temperature at fixed 
filling Ntot/Niat = 0.68, where the green and red curves in- 
dicate the temperature dependence of the superfiuid order <f) 
and checkerboard order 

From the previous discussion, we conclude that onsite 
interactions strongly shift the phase boundary between 
superfluid and supersolid. We now investigate this ef- 
fect in detail at different fillings. We choose a cavity 
detuning Ac ~ —^OOluji, a scattering length of 2.5as and 
a lattice depth Vp = IBEn of the standing- wave pump 
laser, motivated by the recent experiment [8|. Fig. [5] dis- 
plays the resulting checkerboard order 4> (blue line) as a 
function of filling, where four possible phases of the BEC- 
cavity system are observed. Panels (a)-(d) in Fig.[2]show 
the density distribution in real space (left) and in quasi- 
momentum space (right): a) superfluid phase (0 7^ 
and $ = 0) with off-diagonal long-range order (phase 
coherence), b) supersolid ((/> 7^ and $ 7^ 0) with coex- 
isting diagonal long-range order (periodic density mod- 
ulation) and phase coherence, c) Mott insulator (0 = 



and $ = 0) with zero mean-photon number in the cav- 
ity mode, and d) checkerboard solid (</) = and $ 7^ 0) 
with diagonal long-range order and finite mean-photon 
number in the cavity mode. Let us now discuss the un- 
derlying mechanism for the buildup of the self-organized 
phases. The excitation of the cavity mode is a collective 
effect due to all the atoms in the cavity and depends on 
the total particle number, i.e. the more atoms are in the 
cavity, the more photons will be coherently scattered into 
the cavity mode, and the easier the checkerboard pattern 
of the density distribution can be formed. In the absence 
of induced long-range interactions, there are two possi- 
ble phases for strongly interacting bosonic gases in an 
optical lattice: superfluid and Mott insulator. The low- 
lying excitations of the superfluid phase are gapless sound 
modes which can be easily excited while the lowest 
excitations of the Mott insulator are gapped particle-hole 
pairs with an energy gap of order U |24{ . These differ- 
ent excitation properties strongly influence the buildup 
of the self-organized phases. As can be seen from the 
blue curve in Fig. [51 the order parameter $ becomes fi- 
nite with increasing total particle number, and decreases 
to zero again in the vicinity of the Mott insulator. With 
further increase of the filling n > 1, the checkerboard su- 
persolid phase appears again. Interestingly, there is also 
a checkerboard solid phase emerging at n = 1.5, since for 
larger particle number more photons are scattered into 
the cavity mode, and the resulting standing wave in the 
cavity direction suppresses tunneling of atoms and there- 
fore superfluidity. Interestingly, we observe a maximum 
of the order parameter $ at finite temperature due to the 
competition between superfluid and checkerboard order. 
All four phases can be detected experimentally by com- 
bining timc-of-flight measurements and the detection of 
photons leaking from the cavity d, [3l ■ 

We have so far studied the homogeneous case, but in 
real experiments the external trap induces inhomogene- 
ity and a resulting coexistence of superfluid, Mott insu- 
lator, supersolid and checkerboard solid. We will now 
investigate the effect of inhomogeneity on the buildup 
of self-organized phases of the BEC-cavity system, and 
answer the question how the different phases shown in 
Fig. [21 will manifest themselves in the experiment. In 
contrast to the situation with pure contact interactions, 
we find that the properties of the BEG in the optical cav- 
ity are strongly influenced by the trapping potential, due 
to cavity-mediated long-range interactions which are self- 
consistently determined by the density distribution of the 
whole system. Here we consider a A^iat = 32 x 32 lattice 
with harmonic trap strength Vtrap = O-OOSEfi. All other 
parameters are chosen as in Fig. [2l In Fig. [3l we show 
the resulting density (upper panels) and superfluid order 
parameter distributions (lower panels) in real space for 
different total particle numbers. In general, the larger 
the total particle number, the more photons are scat- 
tered into the cavity mode, and thus the easier the sys- 
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buildup of these self-organized phases to be strongly in- 
fluenced by an external trap, due to the density depen- 
dence of scattering between pump laser and cavity mode 
by atoms in the cavity. Self-organized phases can be de- 
tected by combining time-of-flight measurements and the 
detection of photons leaking from the cavity Q , while the 
coexistence of different phases in the presence of an ex- 
ternal trap could be directly observed by quantum gas 



microscopy with single-site resolution [251-127 



FIG. 3. Density distribution n and superfluid order parameter 
versus position on a square (32 x 32) lattice for different 
atom numbers Ntot = 139, 167, 184 and 220 in panels (a), 
(b), (c) and (d), respectively. Other parameters are Ac — 
-500Er, k = 300£;h, Ua = -O.IEr, Vp = 15Er, with a 
harmonic trap Krap = 0.003i5_R. 
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tem can form the self-organized phase. We observe that 
at A^tot = 139, there is almost no checkerboard phase re- 
gion, as visible in panel (a). At Atot = 167, the supersohd 
phase can be clearly observed in the center of the trap, 
since with increasing A^tot the superfluid core expands at 
the trap center and hence more photons are scattered 
into the cavity mode. From Fig. [2j we expect that the 
self-organized phase will disappear again when the num- 
ber of particles increases to a value at which a Mott gap 
arises in the center of the trap, which is clearly visible 
in panel (c) at A'tot = 184. After further increase of the 
particle number to Atot = 220, the checkerboard order 
reappears again. Moreover, we observe a checkerboard 
solid core with average fllling n = 0.5 building up, which 
in the homogeneous case only exists at filling n = 1.5 for 
the same parameters, i.e. we observe the checkerboard 
solid in a parameter regime where it is not visible in the 
homogeneous system. This is because a spatially inhomo- 
geneous atom density due to the harmonic trap strongly 
enhances the scattering between pump laser and cavity 
mode. Observation of these different phases is possible 
by using singe-site addressing techniques in an o ptic al 
lattice based on optical or electron microscopy (25- 27|. 

In conclusion, we have investigated self-organized 
phases (supersolid and checkerboard solid) of both ho- 
mogeneous and trapped ultracold Bose gases coupled to 
a high-finesse optical cavity. We have found that these 
phases are robust against strong onsite interactions at 
zero temperature, where the self-organization phase tran- 
sition is solely driven by quantum fluctuations. We ob- 
serve that thermal fluctuations can enhance the buildup 
of self-organized phases at finite but low temperature. 
In the presence of an external harmonic trap, the co- 
existence of superfluid, Mott-insulating, supersolid and 
checkerboard solid domains is observed. We find the 
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